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Abstract. We use retrodictive quantum theory to analyse two-photon quantum imaging systems. The 
formalism is particularly suitable for calculating conditional probability distributions. 
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1 Introduction 

Two-photon quantum imaging has been studied exten- 
sively for a number of years, both experimentally and the- 
oretically pp. The phenomenon relies upon entanglement 
found in pairs of photons which are produced by sponta- 
neous parametric down-conversion in a x*- 2 ) crystal A 
typical system is shown in figure 1. A pump beam pro- 
duces a pair of entangled photons in a type II downcon- 
version crystal. Due to the properties of both the pump 
field E p and the nonlinear crystal (NLC) these photons 
are entangled in energy and wavevector. The two photons 
within each photon pair are emitted with polarisations 
orthogonal to one another, which enables their separation 
by means of a polarising beam splitter. After reflection 
or transmission at the beam splitter the photons travel on 
their respective paths to be detected at spatially separated 
detection systems. In arm 1 the photon usually propagates 
to a mask of some type (with transmission function t(x)), 
whose image we wish to form, and then after propagation 
it travels to the detector D\ , where it can be recorded at a 
particular position in the transverse plane. In arm 2 there 
is not usually any mask, simply propagation to the detec- 
tor D 2 - It is found that information about the object in 
arm 1 can be found at the detector in arm 2, even though 
the two paths may be widely separated, so that there is 
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no chance that the photon in arm 2 could have interacted 
with the object in arm 1. Of course this can only occur 
when the photon in arm 1 causes the detector D\ to fire, 
so the information is conditional on the occurrence of this 
event. 

A calculation of the spatially-dependent probability 
distribution for joint photodetections at transverse posi- 
tion x\ in arm 1 and position x 2 in arm 2, P(xi,a;2), can 
give information about the object in arm 1. The informa- 
tion is most directly obtained, however, from the condi- 
tional probability distribution that there is a photodetec- 
tion at X2 given that there is one at x\, P{x 2 \x\). This 
probability distribution is what is actually produced at 
D2 in a multi-shot experiment, as the detections in arm 2 
are only recorded if there is also one at detector D\ . It can 
be found from the joint distribution using Bayes' theorem 
0, 



P{x 2 \xi) = , 

P .Tl 



(1) 



Fig. 1. A schematic 2- photon quantum imaging system. Arm 
1 (upper) contains a mask denoted by t between the crystal 
and the detector. This object can be imaged in arm 2. 



where of course, we assume that the arm 1 detection oc- 
curs within a small neighbourhood of x\ in the transverse 
plane, and take the limit that the size of this neighbour- 
hood tends to zero. There is redundant information in the 
joint probability distribution P(x\,X2)- It contains infor- 
mation, for example, about whether any photocounts are 
recorded by either detector due to the fact that the non- 
linear crystal normally does not produce any photon pairs 
within a detector integration time. The conditional prob- 
ability disregards this extra information, as it only deals 
with cases where a photon is recorded at x\. For this rea- 
son it would be better to calculate the conditional prob- 
ability directly but, as we shall see later, there is no way 
to do this in conventional predictive quantum mechanics. 

Klyshko ^ has suggested an advanced-wave interpre- 
tation which has proved useful for the understanding of 
the results of such experiments. In essence the detected 
state in arm 1 is thought of as evolving backwards through 
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Fig. 2. Unfolded version of figure 1. The state evolves from 
the detector on the left towards the detector on the right. 
The Greek letters denote the spatial functions of the 1-photon 
states at each point in the apparatus. 



the system to the crystal, where a conditioning of the 2- 
photon state takes place, forming a 1-photon wavefunction 
which evolves forward in arm 2 and is imaged at the de- 
tector. The situation is similar to figure 2, which shows an 
unfolded version of the system. The state evolves and the 
light is thought of as propagating backwards through the 
system from the detector to the crystal. Then the state 
evolves forward in time as the light propagates from the 
crystal to the detection system in arm 2. 

In this paper we utilise retrodictive quantum theory 
|5lfil7| . that is, quantum theory in which the state of the 
system at any time between preparation and measurement 
is assigned on the basis of a measurement performed on 
the final state rather than the initially prepared state, to 
calculate directly conditional probabilities, such as eq. JTJ, 
for quantum imaging systems. This approach has much 
to recommend it. Only the required probability is calcu- 
lated. The redundant information in the joint distribution 
is not calculated because it is not useful. This is the main 
advantage of the retrodictive approach. Furthermore, it 
provides quantitative predictions based on a formal struc- 
ture in which the reverse-time evolution of the measured 
state corresponds to Klyshko's advanced wave interpreta- 
tion. This supports the Klyshko interpretation, provides 
a formal derivation of conditional probabilities, and thus 
makes the interpretation quantitative. 

The paper is organised as follows. In sectional we de- 
scribe the basic features of retrodictive quantum theory. 
In the following section we apply it to a general quantum 
imaging system such as in figure 1. We then apply the 
theory to a specific example. Finally we summarise our 
results and conclusions. 



where U is the unitary evolution operator which evolves 
the initially-prepared state from the preparation time t p 
to the measurement time t m = t p + r. 

If we do not know which state the preparation device 
prepared, but only have access to the results of the mea- 
surement, then we require not the predictive but the retro- 
dictive conditional probability P(i\j). This is the proba- 
bility that the state pi was prepared given that measure- 
ment result j was recorded. There are two ways in which 
we can calculate this probability. Either we can calculate 
all possible predictive conditional probabilities using pre- 
dictive quantum mechanics, and then use Bayes' theorem 
to find the retrodictive probability, or we can use retrodic- 
tive quantum theory |5lfil7j . Retrodictive quantum theory 
is specifically designed to give the same results as predic- 
tive quantum theory combined with Bayes' theorem 
The Bayesian approach, however, is both more calcula- 
tionally intensive and less elegant than using retrodictive 
quantum theory. 

In retrodictive quantum theory the state of a quantum 
system at any time between preparation and measurement 
is the measured state evolved backwards in time. At the 
preparation time the evolved measured state collapses on 
to the preparation basis. It has been applied to both closed 
systems, in which the time symmetry inherent in quantum 
theory simplifies calculations greatly [5], and to open sys- 
tems, where the retrodictive state evolves backwards in 
time according to a retrodictive master equation analo- 
gous to the Lindblad master equation of predictive quan- 
tum theory [EC]. In closed systems the retrodictive con- 
ditional probability that the prepared state was pi given 
that the later measurement result is j is 



P(i\j) 



Tr (P(i)p lPj (t p )) 
E k Tr(P(k)p k p 3 (t p )) 

Tr( P^PiU^^pjitmpir) 



EfcTr (p(k)p k U^(T)p j (t m )U{T) 



(3) 



2 Retrodictive quantum theory 

Quantum theory is normally formulated in a predictive 
manner. It is particularly useful if we wish to predict 
the outcomes of experiments given particular initially pre- 
pared states. Thus it provides predictive conditional prob- 
abilities. If we have a preparation device which prepares 
states pi with a priori probabilities P(i), and we measure 
these states with a device whose outputs j are describable 
by a probability operator measure (POM) with positive 
elements II j such that Ylj ^ = 1 El j then the predictive 
conditional probability that we obtain the result j if the 
prepared state was pi is given by 

P(j\i) = Tr (pi{t m )fl ) = Tr (#(T)M*p)0V)#i) >( 2 ) 



where the retrodictive state Pj(t m ) — Ilj/Trllj is the nor- 
malised POM element corresponding to the measurement 
result. This evolves backwards in time from the measure- 
ment time to the preparation time, when it collapses on 
to one of the states which could have been prepared. 

It is clear that there is an asymmetry in the forms of 
the predictive and retrodictive conditional probabilities, 
equations J2J an( l ©■ This is not due to any inherent 
time-asymmetry in quantum theory. Rather it is due to a 
choice in standard quantum theory to treat the predictive 
conditional probability as fundamental, and normalise the 
operators which describe prepared and measured states 
differently. Such a choice is not necessary, and when prepa- 
ration and measurement are treated equally the predictive 
and retrodictive conditional probabilities take on symmet- 
ric forms |9I7| . 
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3 Retrodictive analysis 
3.1 General theory 

We now proceed to analyse the general system shown in 
figure 1. We wish to calculate the conditional probability 
distribution of detection at a general transverse position in 
arm 2 given a detection at a particular transverse position 
in arm 1. We will formulate the theory in one transverse 
dimension x. Extension to the whole transverse plane is 
straightforward. 

In conventional quantum theory a fully predictive cal- 
culation is performed based on the two-photon state pro- 
duced by the crystal evolved forward in time and space 
through both paths to form the joint probability distri- 
bution of one detection in each path. We will calculate 
only the conditional probability, performing a calculation 
which is part retrodictive and part predictive in nature, 
in the spirit of the Klyshko interpretation of such exper- 
iments. In order to simplify calculations further we will 
dispense with the formal structure of density operators 
and POM elements describing preparation and measure- 
ments, and simply use prepared and detected states. 

Suppose that a photon is registered by a detector cen- 
tred at transverse position x\ in arm 1. This is represented 
by the 1-photon state 



lb, 



dxa(x)d' (x)\0) 



where a(x) is a normalised complex function of trans- 
verse position centred on x\, so that J dx\a(x)\ 2 = 1. 
This function gives the spatial profile of the detector. The 
continuous-mode annihilation operator a(x) and the con- 
jugate creation operator obey the commutator |1U) 



[a(x), a/{x')\ — 5{x — x'). 



(5) 



The 1-photon retrodictive state can be evolved back- 
wards in time from the detection time. As it does so, we 
can follow the spatial profile back through the apparatus 
to the point of preparation. This approach is typical of 
Fourier optics [TT]. We denote the various functions of x 
at different points in the apparatus by a\,a 2 etc (fig. 2). 
The first part of the propagation is the propagation to the 
object. This is represented by convolution of the spatial 
detector function a(x) with another function of x, h\(x) 
to take account of the propagation. The state is still a 
1-photon state, but its spatial profile has become 



a\{x) — I dx'a{x')h\(x — x'). 



(6) 



The object which is to be imaged is accounted for by a 
simple transfer function t(x), which is a spatially- varying 
complex function whose modulus is not greater than unity. 
Thus 



ct 2 (x) — ai(x)t(x). 



(7) 



Note that a one-photon wavefunction with spatial func- 
tion defined by eq. is not normalised. This is not a 



problem as we simply normalise probabilities at the end 
of the calculation. The next part of the propagation is 
from the object to the crystal. Again this is accounted for 
by convolution 



&3(x) = J dx' a 2 {x')h 2 {x — x') 

fdx'a 1 (x')t(x')h 2 {x - x') 
fdx'[ dx"a{x")h 1 {x' - x")t(x')h 2 (x - x'). (8) 



Thus the retrodictive state at the crystal is the 1-photon 
state of the form defined by eq. (@J , but with spatial profile 
a replaced by the convolution 

We now condition the predictive state of the crystal 
using the retrodictive state from arm 1. The output of the 
crystal is assumed to be a 2-photon state of the form 

|2^)i,2 = fdx f dx'f3(x,x')^(x)b\x')\0) h2 , (9) 



where w is the creation operator for arm 2. On condition- 
ing this forms the one photon state in arm 2 



1 1 a- ) 2 — l(^x\^x',x" )i. 2 = / dx(3 1 (x)P(x)\0} 2 , (10) 



(4) where 



P x (x) = / dx'at{x')P(x',x) 



(11) 



It is clear that by conditioning the 2-photon crystal state 
with the retrodictive state from the detector in arm 1 we 
produce a 1-photon state in arm 2. In fact the combination 
of the detector in arm 1 and the crystal formally constitute 
a quantum state preparation device |9l7j . The complex 
conjugate in this function reflects the fact that the state 
evolution in arm 1 has been backwards in time. 

The final part of the calculation consists of propagation 
to the detector in arm 2. This is again taken account of 
by convolution with the propagation function ha(x). Thus 



(3 2 (x) — I dx'Pi(x')h 3 (x - x 1 ), 



(12) 



with the state given by eq. I|1U|) with (3 2 as the spatial 
profile. 

The conditional detection probability distribution for 
obtaining a detected photon at position x 2 in arm 2 given 
a detection at x\ in arm 1, P(x 2 \x\) is simply given by the 
squared modulus of the final spatial profile, effectively a 
multiple spatial convolution of all of the spatial functions 



P(x 2 \xi) 



\fr{x 2 )Y 



J dx 2 \(3 2 (x 2 )\' 



(13) 



Note that we now must divide by the integral of the func- 
tion in order to normalise the probability distribution. In 
principle we could have renormalised the 1- and 2-photon 
wavefunctions and this would have had the same effect. 
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Fig. 3. The system considered in section 3.1. A pair of photons 
produced by the crystal are separated at a polarising beam 
splitter (b/s) and detected at separate detection systems. 

3.2 Example: Direct imaging of an object 

The utility of this general approach and of the formula 
derived in the previous section can be illustrated by the 
following simple example (see fig. 3). Suppose that there is 
a lens of focal length / placed in each arm of the system. 
The pair of photons produced by the crystal is separated 
using a polarising beam splitter so that only one photon 
can be counted at each detector. The detector in arm 1 is 
placed at the focal length of the lens. The distance from 
the lens to the object to be imaged is arbitrary. We as- 
sume, however, that the object, the beam splitter and the 
crystal are all sufficiently close together that the small 
amount of propagation here has no effect, so h,2(x) = 1. 
In arm 2 the crystal is placed at the focal length of the 
lens, and a spatially-resolving detector is placed after the 
lens. Any propagation after the lens then has no effect on 
the results found at the detector in arm 2. 

The arm 1 detector resolution function will generally 
be of 'top hat' form, but we will use a Gaussian 



"0*0 = ( -^o 

1 1T(T Z 



1/4 



-(x-x l ) 2 /(2a 2 ) 



(14) 



for ease of calculation. In any case when the spatial res- 
olution of the detector is either very good or very poor 
the exact form of the function will not matter. It will be 
useful to write the state of the system in Fourier space 



II, 



1» 



dk x a(k x )a\k x )\0) , 



where 



a(k x ) 



a(k x 



dxa(x)e 



-ik~x 



9\ 1/4 

° X r -k 2 I a 2 /2-ik a: x 1 



dxd(x)e 



-ik x x 



(15) 

(16) 

(17) 
(18) 



are the transverse spatial Fourier transforms of the spatial 
function and operator, and the vacuum state is now the 
state of no photons at any transverse wavevector. For the 
Gaussian detector profile given above the detected state 
has a transverse wavevector profile which is also a Gaus- 
sian. 

Propagation back to the lens corresponds to a modifi- 
cation of the transverse wavevector profile, 



5ii{k x ) 



2\ 1/4 



D ifk z „-Ka 2 /2iklf /k z -ikxX! 



(19) 



The lens effectively takes the Fourier transform of this 
function, so that components which propagate with dif- 
ferent transverse wavevectors between the detector and 
the lens all propagate with the same transverse wavevec- 
tor from the lens to the object, but with spatial profile 
given by 



1 



1/4 



na* J (l-2if/k z a 2 ) 1 / 2 
x exp [-{x - x 1 f/{2a 2 (\ - 2if/k z a 2 ))]. (20) 

As the spatial profile propagates unidirectionally, the dis- 
tance from the lens to the object is arbitrary, and we do 
not consider it. After propagation back through the object 
the spatial profile becomes 0:2(2;) = t{x)a\(x), with ot\(x) 
given by eq. H2U|) . As was stated earlier, we assume that 
the crystal and the polarising beam splitter are placed 
sufficiently close to the object that the small amount of 
propagation involved makes no difference. Then 



otz{x) = 012(2;) = t(x)ai{x). 



(21) 



The spatial profile of the 2-photon state is given by 
the functions 



/3(k X i k x ) 



1 

2^ 



c^[-(k x + k' x Y/2 K % (22) 



0{x,x') = ^6(x-x')e- x ' 2K2/2 . 



(23) 



The spread in transverse wavevector k of this function 
corresponds to the spread in transverse wavevector of the 
Gaussian pump beam. Phase matching then ensures that 
the photon pairs have wavevectors related by eq. I|22l) . 
Projection of the back-propagated retrodictive 1-photon 
state onto this 2-photon state produced by the crystal 
produces a 1-photon state with profile given by eqs. lfl"T|) . 
(JUJl and (T23, 



f3i(x) = y/irt* {x)a\{x)e 



-x 2 k 2 /2 



(24) 



and (3i(k x ) is found by Fourier transformation. This state, 
prepared by conditioning a 2-photon predictive state with 
a single photon retrodictive state, propagates forward in 
arm 2 from the crystal to the lens placed at its focal 
length. This propagation again corresponds to modifica- 
tion of f3i(k x ) to form 



Mkx)=Pi(k x y fk >e lk * f ^ 



(25) 
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The lens does the same as in arm 1, and effectively takes 
the Fourier transform of the function, giving a probability 
distribution which depends upon the transverse coordi- 
nate as in eq. fT3|l . Again, any further propagation from 
the lens to the detector causes no change in the profile. 

The result given in eq. (|25|l can be specialised for par- 
ticular arm 1 detector profiles. In particular, for a narrow 
profile, given by the limit where a — > 0. For a crystal which 
produces a sufficiently broad spread of wavevectors, the 
transverse probability distribution for detection in arm 2 
takes on the form of |t(x)| 2 . Thus the image of an object 
in arm 1 is formed at the detection system in arm 2, even 
though the photon in arm 2 never interacted with arm 1 
at all. 

The other extreme is given by a broad detector in arm 
1. This gives the marginal distribution P(x) which will 
contain no spatial information about the object in arm 1. 
The image is completely washed out by the broad detector. 

Other propagation/detection systems in arm 2 will 
give a different profile. For example if the lens is placed 
a distance 2/ from the crystal, and the detectors are also 
a distance 2/ from the lens then for a 'point' detector in 
arm 1 the probability distribution in arm 2 is the squared 
modulus of the spatial fourier transform of the function 
t(x). 

4 Conclusion 

In this paper we have used retrodictive quantum theory in 
a 2-photon quantum imaging system to calculate the con- 
ditional probability distribution for detection of a photon 
at a particular transverse position in one arm, given a de- 
tection at another particular position in the other arm. 
The retrodictive state evolves backwards in time from the 
detection in one arm to the nonlinear crystal, where it 
conditions the state of the second photon. This condi- 
tioned state then evolves forward in time in the other arm 
and forms the probability distribution. The approach for- 
malises the interpretation of Klysko. We have calculated 
the general probability distribution as a convolution of all 
of the transverse spatial effects in both arm 1 and arm 2, 
and illustrated this with a specific example. 

The advantage of the retrodictive approach over con- 
ventional predictive quantum mechanics is that only the 
required probability distribution is calculated. Much of the 
information in the full predictive probability distribution 
for obtaining two detections at two distinct points in the 
transverse plane, one in each arm is unnecessary. 

In conventional quantum theory the 2-photon state is 
prepared by the crystal, which forms a state preparation 
device. The two detectors, one in each arm are measure- 
ment devices. The part-retrodictive, part-predictive ap- 
proach that we have described here represents the system 
differently. The crystal together with the detector in arm 
1 formally constitute a composite 1-photon state prepara- 
tion device which prepares a photon in arm 2 whose prop- 
erties are determined nonlocally in time both by physical 
processes in the crystal, and the later details of the prop- 
agation to the detector in arm 1. 
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